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ABSTRACT 

The leaf inclination angle distribution (LAD) is an important characteristic of vegetation canopy structure affecting light 
interception within the canopy. However, LADs are difficult and time consuming to measure. To examine possible 
global patterns of LAD and their implications in remote sensing, a model was developed to predict leaf angles within 
canopies. Canopies were simulated using the SAIL radiative transfer model combined with a simple photosynthesis 
model. This model calculated leaf inclination angles for horizontal layers of leaves within the canopy by choosing the 
leaf inclination angle that maximized production over a day in each layer. LADs were calculated for five latitude bands 
for spring and summer solar declinations. Three distinct LAD types emerged: tropical, boreal, and an intermediate tem- 
perate distribution. In tropical LAD, the upper layers have a leaf angle around 35° with the lower layers having horizon- 
tal inclination angles. While the boreal LAD has vertical leaf inclination angles throughout the canopy. The latitude 
bands where each LAD type occurred changed with the seasons. The different LADs affected the fraction of absorbed 
photo synthetically active radiation (fAPAR) and Normalized Difference Vegetation Index (NDVI) with similar rela- 
tionships between fAPAR and leaf area index (LAI), but different relationships between NDVI and LAI for the differ- 
ent LAD types. These differences resulted in significantly different relationships between NDVI and fAPAR for each 
LAD type. Since leaf inclination angles affect light interception, variations in LAD also affect the estimation of leaf 
area based on transmittance of light or lidar returns. 
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1. Introduction 

Vegetation canopies are complex collections of leaves, 
branches, fruits, and flowers. One of the main purposes 
of the canopy is to absorb solar radiation to power pho- 
tosynthesis. The arrangement and quantity of the leaves 
influence how well light penetrates the canopy, affecting 
the distribution of light to the leaves. Important canopy 
characteristics that determine the interception of radia- 
tion include: 

a) leaf area index, 

b) vertical distribution of the foliage, 

c) the leaf inclination angle distribution, 

d) leaf optical properties, 

e) the distribution of the leaves in the canopy, i.e. the 
clumpiness, and 

f) leaf azimuth angle distribution [1]. 


In the organization of the canopy, plants balance light 
interception against other factors, such as water usage 
and temperature control [2,3]. The same canopy charac- 
teristics determining light interception also affect the 
spectral reflectance of the canopy providing a link be- 
tween canopy architecture and remote sensing (e.g. [4, 
5]). 

Leaf orientation with respect to the position of the sun 
is a key factor in determining the amount of light inter- 
cepted by a leaf, and also affects the fraction of incident 
sunlight that penetrates the canopy to lower layers of 
leaves [6,7]. The orientation of a leaf is described by its 
azimuth and inclination angles. The direction from north 
that the leaf points in is the azimuth angle, and the leaf 
inclination angle is a measure of the droop of the leaf 
from horizontal. Often leaf inclination distributions 
within a canopy are described by mathematical functions, 
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such as spherical, erectophile, or extremophile [8,9]. 

Leaf azimuth and inclination angle distribution may 
vary between plant species, growth stage or even time of 
day [8,10-14]. Some plants show heliotropic leaf move- 
ments, with the leaves tracking the sun throughout the 
day. A diaheliotropic leaf, one that keeps its surface per- 
pendicular to the direction of the sun, can intercept over 
50% more solar radiation than a horizontal leaf over a 
day. Conversely, a paraheliotropic leaf, one that keeps 
the leaf surface parallel to the sun, minimizes light inter- 
ception [15,16]. Leaf angles that reduce light interception 
can reduce midday heat loads on the leaf and increase 
water use efficiency [2,3]. Further, since the photosyn- 
thetic light response of leaves saturates at moderate light 
levels [17], leaves that tend to be vertical may maintain 
high photosynthetic levels during periods of high solar 
elevations although intercepting less light than horizontal 
leaves [18]. While at low solar elevations the vertical 
leaves may have an advantage over horizontal leaves by 
intercepting more light when the incident light is not 
saturating [19]. 

Leaf inclination angles also vary with position in the 
canopy. Duncan’s [20] modeling study showed that hav- 
ing erect leaves at the canopy top and horizontal leaves 
in the lower levels of the canopy maximized mid- sum- 
mer mid-latitude plant photosynthesis. Horn’s [18] model 
also predicts that leaves at the top of the canopy will tend 
to be erect, while leaves low in the canopy will tend to- 
ward horizontal to maximize photosynthesis. The general 
results of these modeling studies are supported by obser- 
vations of forest canopies. Measurements of summertime 
leaf inclination angles in several mid-latitude temperate 
forests show the upper canopy mean leaf inclination an- 
gles between 35° and 55° and understory leaf inclination 
angles around 15° [21-24]. 

When viewed from a point on the Earth, the path the 
sun takes through the sky over a day depends on the time 
of the year and latitude of the observer. If canopy leaf 
inclination angle distributions (LAD) are affected by 
solar illumination, then the distribution of inclination 
angles may be related to latitude. Herbert [25] showed 
that mean leaf inclination angles of the arctic rose ( Dryas 
octopetala L.) increased with decreasing latitude over a 
wide range of northern latitudes. He also found leaf in- 
clination angles for subtropical species followed the 
same relationship with latitude as the arctic rose. Bartlett 
et al. [26] studied salt marsh cord grass (Spartina al- 
terniflora L.) over its range in eastern North America 
from Nova Scotia at 46°N south to Florida at 29°N. They 
examined relationships between Spectral Vegetation In- 
dices (SVI), such as the Normalized Difference Vegeta- 
tion Index (NDVI), and biomass and found an abrupt 
discontinuity in the relationships for cord grass that oc- 
curred at approximately 37°30’N. They attribute the 


change to differences in leaf angles; northern canopies 
leaves were more horizontal than the leaves in the south- 
ern canopies. 

The work of Bartlett et al. [26] points out the impor- 
tance of considering variations in leaf inclination angles 
in the interpretation of remotely sensed data. Modeling 
studies have shown the relationship between NDVI and 
the daily fraction of absorbed photosynthetically active 
radiation (fAPAR) to be linear [27] or smooth curves [4, 
5,28]. Studies looking at the effect of different LADs 
found only small variations the NDVI-fAPAR relation- 
ship [4,5]. However, these studies used theoretical dis- 
tribution functions, assumed the leaf angle distribution to 
be the same throughout the canopy, and did not examine 
differences in illumination patterns for varying latitudes 
for different LADs. 

LAD is an important structural feature of plant cano- 
pies, and one of the more adaptable characteristics, yet 
little is known if identifiable LAD patterns occur in na- 
ture, and, if they do, what their spatial and temporal 
variability are. For a given location and time, are there 
optimal LADs? If so, what does the light penetration 
pattern look like in these LADs? Finally, what effect 
might varying LAD patterns have on the remote sensing 
of biophysical variables? This paper tries to make a first 
attempt to address these questions through a simple 
“thought experiment” by joining a canopy radiative 
transfer model with a simple leaf photosynthesis model 
to predict leaf angle distributions that optimize photo- 
synthesis within canopy layers for a variety of locations 
and times. The model also allows an examination of the 
effect of the resulting LADs on the remote sensing of the 
canopy. 

2. Model Description 

The SAIL (Scattering from Arbitrarily Inclined Leaves) 
canopy model is used to simulate the reflectance and 
absorption of radiation by a vegetation canopy [29,30]. 
The SAIL model calculates energy fluxes at any depth in 
the canopy, and can determine a directional reflectance 
for any combination of sun and view angles. The model 
assumes a canopy made up of uniform horizontal layers 
of infinite extent. The characteristics of each layer are 
determined by the amount, orientation, and optical prop- 
erties of the materials occurring in that layer. In the 
model the distribution of the material within each layer is 
assumed to be random. Leaf azimuth angles are also as- 
sumed to be random. In this study the leaf optical proper- 
ties (i.e. leaf spectral reflectance and transmittance) are 
held constant throughout the canopy using baseline val- 
ues from Go ward and Huemmrich [5]. 

For each layer the absorbed photosynthetically active 
radiation (APAR) is calculated from the energy fluxes 
above and below that layer. APAR for a given layer is 
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found using the relation: 

APAR = PARo - PARr - (PARi - PARsr), (1) 

where PARo is the incoming radiation at the top of the 
layer, PARr is the radiation reflected upward from the 
top of the layer, PARi is the radiation transmitted down- 
ward through the layer and PARsr is the radiation re- 
flected upward into the layer from lower layers and soil 
[31]. In the lowest layer of the canopy the term in paren- 
thesis is the radiation absorbed by the soil. The amount 
of instantaneous APAR is determined from the SAIL 
model and is a function of how LAD and leaf area index 
(LAI) for each layer interact with the sun angle. To de- 
termine the daily fAPAR APAR values are summed 
throughout the daylight period and divided by the total 
incoming PAR for the day [5]. 

APAR must be determined for each leaf as an input 
into the photosynthesis model. This is calculated by di- 
viding the total APAR for a layer by the leaf area of that 
layer. Since this produces an average illumination over 
the layer, it is important to have low leaf areas within the 
layer to minimize the effect of shadowing. This study 
uses layers with a LAI of 0.5. This constraint has no ef- 
fect on the physical thickness of the layer since the SAIL 
model has no explicit measurement of depth in it. 

A simple relationship between photosynthesis rate and 
light intensity for a constant carbon dioxide concentra- 
tion is: 

p = l/[tf + 6/(l-Ic)], (2) 

where p is the rate of dry matter production in g-m” 2 -hr _1 , 
I is the intensity of solar radiation in W-nT 2 , and Ic is the 
light compensation point [17]. The parameters a and b 
describe the controlling factors affecting photosynthesis. 
Under intense light, photosynthesis reaches a maximum 
value of Ha. Therefore, a is proportional to the sum of 
the resistances to carbon dioxide diffusion from the air to 


the chloroplasts. The term b/l may be regarded as the 
photochemical resistance, so that b is inversely propor- 
tional to the quantum efficiency in low light levels. In 
this simulation exercise the value of a is taken as 1 
m 2 -hrg _1 , b as 34.9 W-hrg” 1 , and Ic, is 1 W-nT 2 . These 
values remain constant for all layers; that is, there are no 
sun or shade leaves in the model. Incident solar radiation 
is assumed to contain 45% of energy in the photosyn- 
thetically active region of 0.4 to 0.7 micrometers [32]. 
The direct beam radiation reaching the surface of the 
earth is taken to be 1047 W-nT 2 . Incident PAR at the 
canopy top is the product of the total direct beam PAR 
and the cosine of the solar zenith angle for the given lo- 
cation and time of day. The fraction of diffuse incoming 
PAR, if any, is also added to the incident PAR. 

The photosynthesis model uses the APAR per leaf area 
determined from the SAIL model as an input to calculate 
gross production. The instantaneous photosynthesis val- 
ues are integrated to get daily production. Respiration is 
not included in the calculations. The model begins with a 
single layer with an LAI of 0.5 for a given location and 
day of year (see Figure 1). Leaf inclination angles are 
varied between 5° and 85° in 10° increments and the 
daily productivities are compared. The inclination angle 
with the highest production is saved. The next iteration 
consists of two layers, where the top layer has a fixed 
leaf inclination angle of the maximum producer in the 
previous step. The inclination angles of the lower layer 
are varied to find the maximum daily production as be- 
fore. This process is repeated for several layers, stopping 
if the APAR for the leaves in a layer becomes less than 
the light compensation point. There is no attempt in the 
model to maximize total canopy production, but only to 
maximize the daily production of each layer. This as- 
sumption is true if the canopy consists of different plants 
occupying different canopy layers or if there is little or 
no transport of photosynthetic production between 
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Figure 1. A schematic diagram showing the approach used for modeling vertical patterns of leaf inclination angle distribution 
in vegetation canopies. In Step 1, the modeling begins with a single layer in the SAIL model with an LAI of 0.5 over the soil 
background. Leaf inclination angles are varied between 5° and 85° under the range of solar zenith angles that describe the 
daily path of the sun for a given place and date. Daily productivities for each leaf inclination angle are compared and the 
inclination angle with the highest production is saved. In Step 2, a second layer of leaves is added under an upper layer with 
the fixed leaf inclination angle from the previous step (as indicated by the solid box around the layer in the diagram). The 
inclination angles of the lower layer are varied as before to find the maximum daily production. This process is repeated as 
shown in Step 3 and continued for several more layers. 
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branches of an individual plant. 

The SAIL model also allows the calculation of direc- 
tional reflectances from the canopy. Leaf reflectance and 
transmittance were held constant in the simulations. The 
nadir viewing canopy reflectances in visible and near- 
infrared bands are determined for the local sun angle at 
10:00 AM. This approximates the overpass times of 
Terra, Landsat, and morning National Oceanic and At- 
mospheric Administration (NOAA) satellites. The visible 
and near- infrared reflectances are combined in NDVI. 
NDVI is defined as 

NDVI = (/? NIR — P VIS )/ (/Nir + Pv IS ) (3) 

where /^ir is the canopy reflectance in the near-infrared 
band and ycyis is the canopy reflectance in the visible 
band [33,34]. 

3. Results 

The modeled LAD for canopies ranging from the equator 
to 80°N latitude for both the summer solstice and the 
equinox are shown in Tables 1 and 2. In the summer, the 
leaf inclination angles at the top of canopies from the 
equator to 40°N are around 35°. These angles are found 
in the top three layers, below that low light levels result 
in nearly horizontal inclination angles. North of 40°N the 
upper layers have nearly vertical leaves, due to the lower 
solar elevation angles found in high northern latitudes. 

In the spring and fall at the equinox the modeled LAD 
patterns are different than the summer distributions. At 
the equator the leaves are all nearly horizontal through- 
out the canopy, the effect of the canopy maximizing 
production while the sun is directly overhead. At 20 °N 

Table 1. Simulated leaf inclination angles in degrees for 
latitude bands on the summer solstice (solar declination 
23.5°). Each column contains the leaf inclination angles with 
increasing depth in the canopy for a given latitude. The LAI 
column is the cumulative LAI from the top of the canopy. 
The leaf inclination angles vary between 5° for horizontal 
leaves to 85° for vertical leaves. 


LAI 



Latitude 



0° 

20°N 

40°N 

60°N 

80°N 

0.5 

35 

35 

35 

85 

85 

1.0 

35 

45 

35 

85 

85 

1.5 

35 

35 

35 

85 

85 

2.0 

5 

5 

5 

5 

85 

2.5 

5 

5 

5 

5 

85 

3.0 

5 

5 

5 

5 

85 

3.5 

5 

5 

5 

5 

85 

4.0 

5 

5 

5 

5 

85 

4.5 

5 

5 

5 

5 

85 

5.0 

5 

5 

5 

5 

85 


Table 2. Simulated leaf inclination angles in degrees for 
latitude bands on the equinox (solar declination 0°). Each 
column contains the leaf inclination angles with increasing 
depth in the canopy for a given latitude. The LAI column is 
the cumulative LAI from the top of the canopy. The leaf 
inclination angles vary between 5° for horizontal leaves to 
85° for vertical leaves. Table entries containing a dash (-) 
have APAR levels in that layer that were less than the light 
compensation point. 


LAI 



Latitude 



0° 

20 °N 

40°N 

60 °N 

80°N 

0.5 

5 

35 

85 

85 

85 

1.0 

5 

35 

85 

85 

85 

1.5 

5 

35 

85 

85 

85 

2.0 

5 

5 

5 

85 

85 

2.5 

5 

5 

5 

85 

- 

3.0 

5 

5 

5 

85 

- 

3.5 

5 

5 

5 

85 

- 

4.0 

5 

5 

5 

85 

- 

4.5 

5 

5 

5 

- 

- 

5.0 

5 

5 

5 

- 

- 


the LAD pattern is similar to the summer’s LAD. How- 
ever, 40°N the distribution changes from the tropical 
pattern of mostly horizontal leaves to a more boreal pat- 
tern of erect leaves at the top of the canopy. North of 60° 
N the entire canopy consists of vertical leaves. 

A sensitivity analysis was performed to evaluate the 
effects of variations to the initial simulations. Initially, 
the canopy model contained only leaves, and a second 
model run added branches to the canopy. Since branches 
contain non-photo synthetic material, the model was re- 
structured into pairs of layers, where each layer pair had 
an upper layer of leaves and a lower layer of branches. 
Branch area was taken to be a tenth of the leaf area, or 
0.05 per layer. These simulations result in the same leaf 
angle structure as the leaf only simulations described 
above, although the branches intercept between three to 
seven percent of the PAR. 

Changes in solar zenith angle throughout the day are a 
driver for leaf angles in this model, so another variation 
of the simulations involves altering the fraction of inci- 
dent diffuse, i.e. nondirectional, light. If the proportion of 
diffuse light in the incoming radiation is increased, two 
effects are observed in the leaf angle distributions. First, 
the top layers of the tropical canopies become more erect. 
For example, with 60% incoming diffuse light, the top 
layer leaf angle at 40°N in summer changes from 35° to 
55°. The second effect is that leaves begin to become 
horizontal higher in the canopy. Increasing the propor- 
tion of diffuse radiation for the boreal canopies has no 
effect on their leaf angle distribution. 
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The addition of perturbations to the original simula- 
tions in the form of additional canopy components by 
including branches and increasing the fraction of nondi- 
rectional diffuse light were found to have little effect on 
the general patterns of LAD suggesting some robustness 
to the model results. 

The simulated LAD patterns can be evaluated through 
comparison with measured LAD profiles. It is difficult 
and time consuming to measure leaf angles in forests; 
however, some studies have measured leaf angles along 
with depth in the canopy. Hutchison et al. [23] measured 
the architecture of a mature deciduous forest near Oak 
Ridge, Tennessee at approximately 36°N. The forest 
consisted of Quercus alba L., Q. prinus L., Q. velutina 
Lam., Q. falcata Michx., Acer rubrum L. and Lirioden- 
dron tulipifera L. Measurements of leaf inclinations were 
made with a protractor and plumb line. At the top of the 
canopy the leaf inclination was found to be approxima- 
tely 50° and remained near that value within the main 
canopy. In the understory the leaf angle dropped down to 
near 15°. This pattern is similar to the model results, par- 
ticularly if one assumes a fair proportion of the incoming 
light is diffuse in East Tennessee. 

Another study of canopy architecture was made by 
Miller and Lin [22] for a stand of Acer rubrum L. in 
Coventry, Connecticut at approximately 42 °N. They used 
a point drop method; at each sample point a weighted 
line was dropped from above the canopy and every leaf 
touched by the line was measured. They report that the 
leaves at the top of the canopy have a wide distribution 
of leaf angles. In the upper canopy the probability of en- 
countering any given leaf angle between 15° and 75° is 
approximately 0.2. The point drop method may underes- 
timate the fraction of erect leaves. Since erect leaves 
have a small horizontal cross section, they are less likely 
to be intercepted by the vertical sample. The leaf angles 
in the lower canopy were found to be within 15° of hori- 
zontal with a probability of 0.7. 

Ford and Newbould [21] measured leaf angles in a 
coppiced woodland of sweet chestnut ( Castanea sativa 
Mill.) near Kent in the United Kingdom at approximately 
51°N. They found mean leaf inclination angles at the 
canopy top to be about 36°. The inclination angles 
dropped to 15° when the cumulative LAI from the top of 
the canopy was 5. The upper canopy leaf angles match 
the model results for the summer case at 40 °N; however, 
the model indicates that the transition to horizontal 
leaves should occur higher in the canopy than these 
measurements show. 

In the Craigiebum Range, South Island, New Zealand, 
Hollinger [24] measured leaf inclination angles of the 
evergreen mountain beech ( Nothofagus solandri Hook.f.). 
The site was located at 43°08’S. Mean leaf inclination 
angles were 43.3° for the canopy top, 22.1° for the mid- 


canopy, and 17.0° for the lower canopy. The vertical dis- 
tributions of leaf angle for a 46-year-old stand of Japa- 
nese cypress ( Chamaecyparis obtusa Endl.) in central 
Japan decreased exponentially from top to bottom of the 
canopy, ranging in value from greater than 55° to 30.3° 
[35]. Again, these studies show a general pattern similar 
to the simulated summer case at 40 °N. 

These measurements of the vertical patterns leaf an- 
gles within forest canopies do not verify the results of the 
simulations. However, they show similar patterns to the 
simulated LAD and show leaf angles varying with depth 
in the canopy. 

Using the leaf angle distributions predicted for each 
latitude, the SAIL model calculates the canopy reflec- 
tance in the visible and near-infrared bands and these 
values are used to determine NDVI. The model also cal- 
culates total daily fAPAR as it varies with increasing 
LAI, allowing an analysis of the relationship between 
canopy structure, fAPAR, and NDVI as they vary over 
space and time. 

Figures 2(a) and (b) show how daily fAPAR varies 
with LAI for five latitudes on the summer solstice (solar 
declination 23.5°) and the equinox (solar declination 0°). 
LAI was varied by adding layers to the canopy with the 
leaf angles predicted as described above, so this plot can 
also be interpreted as a measure of cumulative fAPAR 
with depth in the canopy. In all cases, the fAPAR curves 
vary smoothly with LAI. The curves for 0° through 60°N 
on the summer solstice and 0° through 40 °N on the equi- 
nox are nearly identical and saturate at about the same 
value, approximately 97%. The effect of varying solar 
zenith angles appears at higher latitudes. The curve for 
80°N on the summer solstice is the same as the curve for 
60°N on the equinox. The very low sun angles at 80°N on 
the equinox result in a very high percentage of fAPAR 
even with low LAI. 

The relationships between LAI and NDVI at 10:00 
A.M. solar time (Figures 3(a) and (b)) are more variable 
than the relationships between LAI and fAPAR. There 
are different curvatures for most lines, the saturation 
points vary for most curves and there are discontinuities 
in some curves. In the plots for the summer solstice 
(Figure 3(a)) the latitudes between 0°N and 40 °N have 
canopies with the tropical leaf angle distributions and 
show identical relationships between NDVI and LAI. At 
60 °N the relationship shows a discontinuity when the 
canopy reaches a layer where the more erect leaves 
change to horizontal leaves at a LAI of 2. The line for the 
80°N boreal canopy has its own track, having much 
lower NDVI values than the tropical canopies at low LAI. 
On the equinox (Figure 3(b)), the relationships for every 
latitude has its own distinct curve. The curves for 0°N 
and 20 °N are nearly the same as on the solstice but they 
are no longer identical. The curve for 40 °N is nearly the 
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Figure 2. (a) Relationship between canopy LAI and daily 
fAPAR for canopies predicted by the model at five northern 
latitudes on the summer solstice (solar declination 23.5°); (b) 
Relationship between canopy LAI and daily fAPAR for 
canopies predicted by the model at five northern latitudes 
on the equinox (solar declination 0°). 

same as the 60 °N curve on the solstice, showing the same 
discontinuity at a LAI of 2. The 60 °N curve for the equi- 
nox is similar to the 80°N curve on the solstice. The 80°N 
curve for the equinox has the highest values of NDVI. 

The relationship between NDVI and fAPAR produces 
three distinct curves depending on canopy type (Figures 
4(a) and 4(b)). The tropical canopies, 0° through 40 °N on 
the solstice and 0°N and 20°N on the equinox, have an 
upward curvature. While the boreal canopies, 80°N on 
the solstice and 60°N and 80°N on the equinox, have a 
downward curvature. The canopies at 60 °N on the sol- 
stice and 40 °N on the equinox begin with the boreal 
curves and make a transition to the tropical curves at 
higher values of canopy LAI. Figures 4(a) and (b) show 
how attempting to predict percent fAPAR from NDVI 
may result in errors of as much as 40%, if one does 


Figure 3. (a) Relationship between canopy LAI and nadir 
view NDVI at 10 AM solar time for canopies predicted by 
the model at five northern latitudes on the summer solstice 
(solar declination 23.5°); (b) Relationship between canopy 
LAI and nadir view NDVI at 10 AM solar time for canopies 
predicted by the model at five northern latitudes on the 
equinox (solar declination 0°). 

not take into account the effects of canopy structure and 
latitude. 

The variation of leaf angle with depth in the canopy 
may also affect other approaches to remote sensing of 
vegetation. For example, lidars can directly measure the 
vertical distributions of intercepted surfaces in a canopy, 
and from that infer the distribution of leaf area [36]. The 
conversion between the area intercepted by the vertical 
beam and leaf area depends on leaf inclination angles. 
Figure 5 shows the expected interception of a vertical 
beam by layers of a canopy for the 40 °N case. In this 
example, each layer has the same LAI, but there is a sig- 
nificant difference in the interception with depth in the 
canopy due to variations in leaf inclination angles alone. 
Seasonal and spatial variations in leaf angles as shown in 
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(b) 


Figure 4. (a) Relationship between nadir view NDVI at 10 
AM solar time and daily fAPAR for canopies predicted by 
the model at five northern latitudes on the summer solstice 
(solar declination 23.5°); (b) Relationship between nadir 
view NDVI at 10 AM solar time and daily fAPAR for cano- 
pies predicted by the model at five northern latitudes on the 
equinox (solar declination 0°). 

these simulations can result in significant differences in 
the interception of the layers that must be accounted for 
when processing lidar data to accurately describe the 
vertical distribution of leaf area. 

4. Conclusions 



Figure 5. The fraction of a vertical beam intercepted by 
each canopy layer for the 40°N case on the solstice with 
solar declination 23.5° (black bars), and the equinox with 0° 
declination (gray bars). For the canopy layers, 1 is the top- 
most layer, and each layer has an LAI of 0.5. 


tudes, the optimal LAD changes with season, blending 
these two patterns. There are too few measurements of 
LAD, particularly of tropical and boreal canopies, to al- 
low a verification of the model results. 

The model results further indicate that spatial and 
temporal variation of LAD complicates the interpretation 
of remotely sensed data. The simulations indicate rela- 
tionships between NDVI and fAPAR are nonlinear with 
curvature changes depending on LAD and location. The 
tropical and boreal LADs each produce a distinct NDVI - 
fAPAR relationship. Significant errors in estimation of 
fAPAR may occur by not taking into account spatial 
variability in LAD. In the middle latitudes errors in fA- 
PAR estimation from NDVI may also occur if seasonal 
changes in LAD are not taken into account. LAD pat- 
terns are an important, but little understood, factor in the 
analysis of global remote sensing observations. 

This simple modeling “thought experiment” produced 
interesting results, pointing to the need for further work 
using models with more complex 3 -dimensional descrip- 
tions of vegetation canopies along with measurements of 
leaf angles with depth in canopies. 


This simulation exercise suggests that there may be 
global patterns of LADs in vegetation canopies. For 
canopies optimizing photosynthesis within each layer, 
the LADs are driven by location and sun angle resulting 
in latitudinal and seasonal patterns of canopy structure. 
The optimal LAD for tropical regions has leaves in the 
upper canopy that are around 35° and horizontal in the 
lower canopy. In the high latitudes, the optimal canopy 
LAD consists entirely of vertical leaves. In the mid-lati- 
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